Possible proximity of the Mott insulating Iridate Na 2 Ir0 3 to a topological phase: 
Phase diagram of the Heisenberg-Kitaev model in a magnetic field 
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Motivated by the recent experimental observation of a Mott insulating state for the layered Iridate Na2lrC>3, 
we discuss possible ordering states of the effective Iridium moments in the presence of strong spin-orbit coupling 
and a magnetic field. For a field pointing in the (111) direction - perpendicular to the hexagonal lattice formed 
by the Iridium moments - we find that a combination of Heisenberg and Kitaev exchange interactions gives rise 
to a rich phase diagram with both symmetry breaking magnetically ordered phases as well as a topologically 
ordered phase that is stable over a small range of coupling parameters. Our numerical simulations further 
indicate two exotic multicritical points at the boundaries between these ordered phases. 
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In the realm of condensed matter physics, spin-orbit cou- 
pling has long been considered a residual, relativistic cor- 
rection of minor relevance to the macroscopic properties of 
a material. In recent years this perspective has dramatically 
changed, especially due to the theoretical prediction and sub- 
sequent experimental observation of fundamentally new states 
of quantum matter, so-called topological insulators JT], that 
are solely due to the effect of spin-orbit coupling. The topo- 
logical insulators experimentally realized so far are semicon- 
ductors, whose physical properties can be largely captured by 
band theory of non-interacting electrons. It is an interesting 
challenge, for both theory and experiment, to identify an even 
broader class of materials where this physics plays out even in 
the presence of interactions and strong correlations |2|. Good 
candidate materials for the latter are the Iridates EH). These 
5d transition metal oxides are prone to exhibit electronic cor- 
relations and form (weak) Mott insulators, while the relatively 
large mass of the Iridium ions ( Z = 77) gives rise to a com- 
parably strong spin-orbit coupling, which has been found to 
be as large as A « 400 meV [5]. The most common va- 
lence of the Iridium ions in these materials is Ir 4+ . The d- 
orbitals of this 5d 5 configuration are typically split by the sur- 
rounding crystal field, and for the octahedral geometry of the 
IrOe oxygen cage, result in an orbital configuration where five 
electrons occupy the lowered, threefold degenerate t2 g level. 
Spin-orbit coupling will further lift this degeneracy of the t 2g 
orbitals and for strong coupling the effective I = 1 orbital 
angular momentum [Q is combined with the s = 1/2 spin 
degree of freedom carried by the hole of this partially filled 
*2g orbital configuration. This leaves us with two Kramers 
doublets of total angular momentum j = 3/2 and j = 1/2, 
of which the former is of lower energy and fully occupied by 
four electrons, while the partial filling of the latter gives rise 
to an effective spin- 1/2 degree of freedom. 

In this manuscript we focus on the Iridate Na2lr03, in 
which Nali^Oe slabs are stacked along the crystallographic 
c-axis, and the Ir 4+ ions in the layers form a hexagonal lat- 
tice (4). Recent measurements of the magnetic susceptibil- 



ity provide evidence of effective spin- 1/2 moments and mag- 
netic correlations below Tn ~ 15 K indicating that Na2lr03 
is indeed a Mott insulator [4|. Theoretically, it has been ar- 
gued 00 that the interactions between the effective Iridium 
moments in the Mott regime are captured by a combination 
of isotropic and highly anisotropic exchanges, which can be 
tracked back to the spin and orbital components of the effec- 
tive momenta. A microscopic Hamiltonian interpolating be- 
tween these two types of exchanges is given by 

7 — links 

where the <jj denote the effective spin- 1/2 moment of the Ir 4+ 
ions, 7 = x,y,z indicates the three different links of the 
hexagonal lattice, and < a < 1 parametrizes the relative 
coupling strength of the isotropic and anisotropic exchange 
between the moments. For a = the Hamiltonian reduces 
to the ordinary Heisenberg model, while in the opposite limit 
of highly anisotropic exchanges (a = 1) the system corre- 
sponds to the Kitaev model |9|. The latter is known to exhibit 
a gapless spin-liquid ground state (for equal coupling along 
the links) that can be gapped out into a topological phase with 
non-Abelian quasiparticle excitations by certain time-reversal 
symmetry breaking perturbations J9). One such perturbation 
is a magnetic field pointing in the (111) direction, perpendic- 
ular to the honeycomb layer 

-HtlK+h = -HhK — ' ®i ■ (2) 

i 

The main result of our manuscript is the rich phase diagram 
of this model, shown in Fig. [T] Besides two conventional, 
magnetically ordered phases we find a topologically ordered 
phase and two multicritical points, which we will discuss in 
detail in the remainder of the manuscript. 

Numerical simulations- We determine the ground-state 
phase diagram of Hamiltonian |2]) by extensive 'quasi-2D' 
density-matrix renormalization group (DMRG) ifTUl calcula- 
tions on systems with up to N — 64 sites. In particular, 
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FIG. 1: (color online) Ground-state phase diagram of the 
Heisenberg-Kitaev model |TJ in a (111) magnetic field of strength 
h. Interpolating from the Heisenberg (a = 0) to Kitaev (a = 1) 
limit for small field strength, a sequence of three ordered phases is 
observed: a canted Neel state for a < 0.4, a canted stripy Neel state 
illustrated in Fig.[2j;) for 0.4 < a < 0.8, and a topologically ordered 
state for non-vanishing field around the Kitaev limit. All ordered 
phases are destroyed for sufficiently large magnetic field giving way 
to a polarized state. 
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FIG. 2: (color online) a) The honeycomb lattice spanned by unit 
vectors Si — (1, 0) and 3,2 = (1/2, \/3/2). Illustration of magnetic 
states with b) Neel order and c) stripy Neel order. 



we consider clusters of size N — 2 x N\ x N2, which are 
spanned by multiples Nidi and N2S2 of the unit cell vectors 
ai = (1,0) and ti 2 = (1/2,^3/2) as illustrated in Fig. [2] 
It should be noted that the numerical analysis of Hamiltonian 
(j2j is a challenging endeavor, since not only the entire Hilbert 
space needs to be considered (due to the lack of SU(2) invari- 
ance), but one also has to work with complex data types (due 
to the (111) orientation of the magnetic field). Our DMRG 
calculations keep up to m = 2048 states, which is found 
to give excellent convergence with typical truncation errors 
of less than 10~ 8 . We further use periodic boundary condi- 
tions in both lattice directions, which reduces finite-size ef- 
fects. We have determined the phase boundaries in Fig. |T|by 
extensive scans of the ground-state energy, magnetization, and 
their derivatives in the {a, /i)-parameter space ifTTI . 

Magnetically ordered states.- We start our discussion of 
the phase diagram shown in Fig.[TJby first recapitulating previ- 
ous results [8] for the Heisenberg-Kitaev model ([TJ) in the ab- 



sence of a magnetic field. Interpolating the relative coupling 
strength a between the isotropic Heisenberg limit (a = 0) 
and the highly anisotropic Kitaev limit (a = 1) a sequence 
of three phases has been observed (8): The Neel ordered state 
of the Heisenberg limit is stable for a < 0.4, when it gives 
way to a 'stripy' Neel ordered state illustrated in Fig.|2]which 
covers the coupling regime 0.4 < a < 0.8. In the extended 
parameter regime 0.8 < a < 1 the collective ground state is 
a gapless spin liquid. Near a = 1, perturbation theory reveals 
that the gapless excitations of this phase are emergent Majo- 
rana fermions forming two Dirac cones in momentum space. 

Including a magnetic field in the (111) direction a rich 
phase diagram evolves out of this sequence of three phases. 
For the magnetically ordered states we find that the orientation 
of the order in the Neel and stripy AFM phase cants along the 
(111) direction. To further characterize these canted states, 
it is helpful to analyze the independent symmetries of Hamil- 
tonian |2}. Besides the lattice translational symmetry T and 
a reflection symmetry / around the centers of the hexagons, 
there is an additional C3 symmetry, which is a combination 
of a three-fold rotation around an arbitrary lattice site and a 
three-fold spin rotation along the (111) spin axis [12|. Both 
canted phases break a subset of these discrete symmetries of 
the Hamiltonian. The canted Neel order breaks the C| and 
the / symmetries, which thus leads to a six-fold ground-state 
degeneracy in this phase. The canted stripy phase breaks both 
the C3 and translational symmetry (since the ordering pattern 
doubles the unit cell). As a consequence, we also find a six- 
fold ground-state degeneracy in this phase. 

For sufficiently large magnetic field, the order of both 
canted phases is destroyed and they give way to a simple po- 
larized state. Our numerical simulations strongly suggest that 
the transitions between the polarized state and these canted 
states are continuous, which is in agreement with their spon- 
taneous symmetry breaking. On the other hand, the transition 
between the two canted states at finite field strength (indicated 
by the bold line in Fig. |TJ> is found to be first-order. In our sim- 
ulations this is indicated by a sharp drop of the first derivative 
of the energy dE/da as a function of the coupling parame- 
ter a across this transition - as shown in Fig.|3]for increasing 
strength of the magnetic field h. Approaching the endpoint 
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FIG. 3: (color online) Energy jump along the first-order transition 
between the canted Neel and stripy AFM. 
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of this first-order line around h c ~ 0.7 this drop smoothly 
vanishes, which indicates that this endpoint possibly is a tri- 
critical point at which the two canted magnetically ordered 
phases and the polarized phase meet. The existence of such 
a tricritical point can be understood within a Landau descrip- 
tion with two distinct order parameters - corresponding to the 
discrete symmetry breaking of the two different magnetically 
ordered phases - opening a gap to magnon excitations. 

Topological phase- We now turn to the spin liquid phase 
found for coupling parameters 0.8 < a < 1. For the Kitaev 
limit (a = 1) it has previously been argued [9| that an in- 
finitesimal field along the (111) direction will drive the system 
from the gapless spin liquid into a gapped non-Abelian topo- 
logically ordered phase. As we will discuss in the following 
our numerical simulations allow us to confirm the existence 
of such a topologically ordered state for small magnetic field 
strengths not only in the Kitaev limit, but for the full extent 
of the gapless spin liquid phase, as indicated in the phase di- 
agram of Fig. [T] We will further present an independent and 
non-perturbative way to determine the topological nature of 
this phase. This complements the original argument by Ki- 
taev (5i|, which was primarily based on a perturbation expan- 
sion showing that the leading order effect of a small magnetic 
field h is to introduce a topological mass term for the Majo- 
rana fermions - however, such a perturbative argument should 
be carefully tested when applied to a gapless state. For this 
purpose, we consider an additional three-spin exchange term 
k, indicated by the blue bonds in Fig. Eh), in our Hamiltonian 



ijk 



(3) 



In the Kitaev limit (a = 1, h = 0) this Hamiltonian is exactly 
solvable in the same Majorana fermion representation used in 
the solution of the unperturbed Kitaev model (5J. In particu- 
lar, one can prove that the three-spin exchange n breaks time- 
reversal symmetry and gaps out the spin liquid phase into a 
topologically ordered state with non-Abelian excitations, so- 
called Ising anyons. To demonstrate that a small magnetic 
field in the (111) direction drives the system into the same 
phase, we have numerically calculated the phase diagram in 
the presence of both perturbations as shown in Fig. |4] The 
phase boundaries were again obtained by scanning the deriva- 
tives of ground-state energy and magnetization in the (h, k)- 
parameter space. In particular, this phase diagram shows that 
one can adiabatically connect the phase for large n and van- 
ishing magnetic field with the phase for small, non-vanishing 
magnetic field and n = 0, thus proving that the magnetic field 
gaps out the spin liquid into the same non-Abelian topological 
phase stabilized by the three-spin exchange. The only feature 
in the diagram is a single phase transition line which sepa- 
rates the topologically ordered state from the fully polarized 
state expected for large magnetic field strengths. For the Ki- 
taev limit (k = 0) this transition occurs for h c ~ 0.072. It is 
interesting to note that the critical field h c initially grows with 
increasing k, but then saturates to some finite value around 
k > 6. Physically, this saturation can be understood by the 
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FIG. 4: (color online) Top panel: Ground-state phase diagram of the 
Kitaev model (a = 1) in the h — k plane, where h is the strength of 
a magnetic field pointing in the (111) direction and k is the strength 
of a time-reversal symmetry breaking three-site term. Lower panel: 
Magnetization sweeps and its derivative for various k. 



behavior of the gap for the Majorana fermions in the exact 
solution for h = 0. The dispersion of the Majorana fermion 



is given by E^ = 2 



1 



+ k 2 sin 2 (&; • Si). 



For small k<1, the Majorana fermion has a gap E g 



However, for large k S> 1 the gap of Majorana fermion re- 
mains finite and independent from k, given by E g ~ 2. Since 
the magnetic field strength h c required to destroy the topo- 
logical phase is determined by the Majorana fermion gap at 
h = 0, the critical field h c thus also increases and then satu- 
rates at large k. 

Field-driven transition out of the topological phase.- We 
now return to the phase diagram of the Heisenberg-Kitaev 
model in Fig. [T] and focus on the transition between the topo- 
logically ordered state and the polarized state for large field 
strength. For the Kitaev limit {a = 1) this transition occurs 
at a critical field strength of h c ps 0.072 and remains almost 
constant as the coupling parameter a is decreased. Interest- 
ingly, our numerics suggest that this field-driven phase transi- 
tion might be continuous or weakly first-order. In particular, 
we find that the second-derivative of the ground state energy 
—d 2 E/dh 2 at this transition diverges with increasing system 
size, while the magnetization M(h) does not show any dis- 
continuity, as shown in Fig. [5^) and b), respectively. 

While the limited system sizes in our study do not allow 
to unambiguously determine the continuous nature of this 
field-driven phase transition, our numerics nevertheless pro- 
vide some further insights what might cause such a continuous 
transition lfl3l . To this end, we plot the number of vortices in 
the ground state as a function of magnetic field, i.e. the num- 
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. 5: (color online) Phase transition for the Kitaev model (a — 1) 
in a (111) magnetic field of strength h. a) Second derivative of the 
ground state energy —d 2 E/dh 2 , b) Magnetization M(h) and its first 
derivative dM(h) /dh for different system sizes, c) Vortex number. 



ber of plaquettes with a non-trivial flux, in Fig. |5j;). Below 
the critical magnetic field, i.e. h < h c , there are no vortices 
indicating a deconfined phase as expected in the presence of 
a vortex gap. At the phase transition, however, the vortices 
appear to condense and the number of vortices in the ground 
state quickly increase above the critical field strength. The 
nature of the phase transition might thus be framed in terms 
of a confinement-deconfinement transition of a non-Abelian 
gauge field, akin to the confinement-deconfinement transition 
in the Abelian discrete gauge theory T\4\ . An example of the 
latter is the Z2 gauge theory, e.g. the toric code in a magnetic 
field Q3) , for which it is well known that flux condensation 
leads to a confinement transition |[T6l . 

A second multicritical point- Finally, we note that there 
appears to be a second multicritical point in our phase diagram 
around a ps 0.8 and h = 0, where the stripy AFM phase and 
the gapless spin liquid meet. We find that in the presence of 
the magnetic field the transition lines of the field-driven phase 
transition out of the corresponding canted and topologically 
ordered states bend in and merge only in the zero-field limit 
as depicted in the inset of Fig. [T] To show that there is in- 
deed no direct transition between the canted stripy Neel state 
and the topological phase we have made extensive scans in 
the coupling parameter a in the vicinity of this putative mul- 
ticritical point for small field strength. As shown in Fig.|6]for 
h = 0.06, the second derivative d 2 E/da 2 of the ground-state 
energy clearly shows two peaks proliferating with increasing 
system size indicative of two well separated phase transitions. 
However, the underlying effective theory for such a multicrit- 
ical point is not known, and will be left for further study. 
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FIG. 6: (color online) Constant field scan in the vicinity of the (multi- 
critical) point separating the stripy AFM from the topological phase. 



Outlook- Having established the rich phase diagram of 
the Heisenberg-Kitaev model in a magnetic field, it is interest- 
ing to speculate where one would place the Iridate Na2irC»3. 
While experiments [4 1 report indications of an AFM ordered 
ground state below Tn ~ 15 K, the precise nature of the or- 
der remains open. Given the considerable suppression of the 
ordering temperature Tn in comparison with the Curie-Weiss 
temperature Ocw ~ 116 K [4|, which is typically interpreted 
as an indicator of frustration, an alternative explanation would 
be the proximity to a quantum critical point, such as the mul- 
ticritical point a ps 0.8 in the context of our phase diagram. 
This would bring the material in close proximity to the spin 
liquid phase for a > 0.8 and the topological phase found for 
a magnetic field pointing in the (111) direction. To further 
substantiate this possibility, it is desirable to study the finite- 
temperature phase diagram of our model system and to con- 
sider the effects of disorder, such as site mixing between the 
Ir and Na sites [4 |. Finally, it would be interesting to bring the 
Mott physics discussed in this manuscript in competition with 
the topological insulator phase suggested in Ref . [17] 
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